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Abstract 

At sufficiently large proton energies, Glauber multiple-scattering theory offers good 
opportunities for describing the final state interactions in electro-induced proton 
emission off nuclear targets. A fully unfactorized relativistic formulation of Glauber 
multiple-scattering theory is presented. The effect of truncating the Glauber multiple- 
scattering series is discussed. Relativistic effects in the description of the final-state 
interactions are found not to exceed the few percent level. Also the frequently 
adopted approximation of replacing the wave functions for the individual scattering 
nucleons by some average density, is observed to have a minor impact on the results. 
Predictions for the separated 4 He(e, e'p) response functions are given in quasi-elastic 
kinematics and a domain corresponding with 1 < Q 2 < 2 (GeV) 2 . 
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1 Introduction 

From early in the seventies till recent years, a systematic survey with the 
exclusive A(e, e'p) reaction on a whole range of target nuclei revealed that 
the momentum distributions of bound protons in nuclei are in line with the 
predictions of the nuclear mean-field model. The occupation probabilities for 
the single-particle levels, on the other hand, turned out to be substantially 
smaller than what could be expected within the context of a naive mean- 
field model [1]. This observation provided sound evidence for the importance 
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of short- and long-range correlations for the properties of nuclei [2]. From 
1990 onwards, the scope of A(e, e'p) reactions has been widened. These days, 
rather than for studying mean-field properties, electro-induced single-proton 
knockout off nuclei is used as an experimental tool to learn for example about 
relativistic effects in nuclei [3] or to study fundamental issues like possible 
medium modifications of protons and neutrons when they are embedded in a 
dense hadronic medium like the nucleus [4,5]. Another issue is the question at 
what distance scales quark and gluons become relevant degrees of freedom for 
understanding the behavior of nuclei. Here, searches for the onset of the color 
transparency phenomenon in A(e, e'p) reactions play a pivotal role [6]. In these 
studies one looks for departures from predictions for proton transparencies 
from models using standard nuclear-physics wave functions combined with 
the best available tools for describing final-state interactions (FSI). 

Traditionally, the results of exclusive A(e, e'p) data-taking have been inter- 
preted with the aid of calculations within the framework of the non-relativistic 
distorted wave impulse approximation (DWIA) [7,8]. In such an approach, 
the electromagnetic interaction of the virtual photon with the target nucleus 
is assumed to occur through the individual nucleons, an assumption which 
is known as the impulse approximation (IA). In the most simple DWIA ver- 
sions, an independent particle model (IPM) picture is adopted and the initial 
and final A-nucleon wave functions are taken to be Slater determinants. The 
latter are composed of single-particle wave functions which are solutions to 
a one-body Schrodinger equation. Typically, in a DWIA approach the final 
proton scattering state is computed as an eigenfunction of an optical poten- 
tial, containing a real and imaginary part. Parameterizations for these optical 
potentials are usually not gained from basic grounds, but require empirical 
input from elastic p + A — > p + A measurements. In the DWIA calculations, 
one adopts the philosophy that potentials which parameterize FSI effects in 
elastic p + A — > p + A reactions, can also be applied to model the ejectile's 
distortions in A + e — > A — 1 + e' + p. 

A concerted research effort which started back in the late eighties has re- 
sulted in the development of a number of relativistic DWIA (RDWIA) models 
for computing A(e,e'p) observables [9,10,11,12,13]. These theoretical efforts 
very much followed the trend of developing relativistic models for p + A — > 
p' + A processes as a potential improvement to the traditional non-relativistic 
distorted- wave approaches. Along the lines of the DWIA approaches, the RD- 
WIA frameworks are usually formulated within the context of the independent- 
particle approximation with wave functions of the Slater determinant form. 
Relativistic bound-state single-particle wave functions are customarily ob- 
tained within the framework of the Hartree approximation to the a — u model 
[14] . Scattering states by solving a time-independent Dirac equation with rela- 
tivistic optical potentials. Systematic surveys illustrated that in some A(e, e'p) 
observables relativistic effects are sizeable [15,16] and that RDWIA-based cal- 
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culations of A(e, e'p) observables are at least as successful as the more tradi- 
tional non-relativistic DWIA theories. 



For proton lab momenta exceeding roughly 1 GeV, the use of optical poten- 
tials for modeling FSI processes does not seem to be natural in view of the 
highly inelastic character and diffractive nature of the underlying elementary 
nucleon-nucleon scattering cross sections. In this energy regime, an alternative 
description of FSI processes is provided in terms of Glauber multiple-scattering 
theory [17,18,19]. This theory essentially relies on the eikonal approximation 
and the assumption of consecutive cumulative scattering of a fast proton on a 
composite target containing A-l "frozen" point scatterers. Most A(e, e'p) cal- 
culations based on the Glauber approach are non-relativistic and factorized. A 
non-relativistic Glauber model for modeling FSI effects in exclusive 4 He(e, e'p) 
reactions was recently pointed out in Ref. [20] . A non-relativistic study of the 
Glauber formalism in d(e, e'p)n can for example be found in Refs. [21] and 
[22]. For (e, e'p) reactions off nuclei heavier than 12 C and nuclear matter, non- 
relativistic Glauber calculations have been reported in Refs. [23,24,25,26]. 



In this work we wish to present a relativistic formulation of Glauber theory for 
calculating A(e, e'p) observables. The major assumptions underlying our rela- 
tivistic and unfactorized model bear a strong resemblance with those adopted 
in the RDWIA models developed during the last two decades. One of the 
pivotal assumptions underlying Glauber theory is the eikonal approximation. 
Expressions of Dirac-eikonal scattering amplitudes for elastic p + A processes 
have been derived in Ref. [27] and for A(e,e'p) reactions in Refs. [28,29]. 



The organization of this paper is as follows. First, we briefly sketch the A(e, e'p) 
formalism in Sect. 2. Next, we introduce the relativistic formulation of Glauber 
multiple-scattering theory in Sect. 3. Section 4 is devoted to a presentation of 
results for the Dirac-Glauber phase for the nuclei 4 He, 12 C, 56 Fe and 208 Pb. 
The Dirac-Glauber phase is a function which accounts for all FSI effects when 
computing the A(e, e'p) observables. The contribution of single- and multiple- 
scatterings is estimated for the target nuclei 4 He, 12 C and 208 Pb. Attention is 
paid to the role of relativistic effects when computing the impact of final-state 
interactions. The validity of a frequently adopted approximation, namely the 
replacement of the individual nucleon wave functions by some average nuclear 
density, is investigated. In Sect. 5 we present predictions of our model for 
the separated (e, e'p) response functions for a 4 He target nucleus. Section 6 
summarizes our findings and states our conclusions. 
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2 A(e, e'p) Formalism 



We follow the conventions for the A(e, e'p )B kinematics and observables in- 
troduced by Donnelly and Raskin in Refs. [30,31]. The four- moment a of 
the incident and scattered electron are denoted as Kg(e,k) and K£(e',k'). 
The electron momenta k and k' define the lepton scattering plane. The four- 
momentum transfer is given by (u, q ) = = K£ — K£ = K\_ x + K£ — K%, 

where K^(Ea,]ca), ^a-^C^a-i, &a-i) and K£(E p ,k p ) are the four-momenta 
of the target nucleus, residual nucleus and the ejected proton. The z-axis lies 
along the momentum transfer q and the y-axis lies along kxk'. The hadron re- 
action plane is defined by k p and q. The electron charge is denoted by — e, and 
we adopt the standard convention Q 2 = —q^ > for the four-momentum 
transfer. 



In the laboratory frame, the differential cross section for A + e — > (A — 1) + 
e' + p processes can be computed from [30,31,32] 
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where r] e (K e >, S"; K e , S)^ is the electron tensor, and the nuclear tensor W^ v 
is given by 
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Here, J M is the electromagnetic current operator, | A(0 + , g.s.)) the ground state 
of the target even-even nucleus and \ A — 1(J R Mr)) the discrete state in which 
the residual nucleus is left. 



For m e /e <C 1 the contraction of the electron tensor 77^ with the nuclear one 
W^ u results in an expression of the form [30] 

4m 2 e r) e (K e ,, S'; K e , S) pv W^ = v E vkKk , (6) 
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where the label K takes on the values L, T, TT and TL and refers to the 
longitudinal and transverse components of the virtual photon polarization. 
The TZk are the nuclear response functions and contain the nuclear structure 
and dynamics information. Further, vq = (e + e') 2 — q 2 and the Vk depend 
on the electron kinematics. Combination of the above results leads to the 
following final expression for the A(e, e'p) differential cross section [30,31] 
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where % is the Mott cross section 
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8 e the angle between the incident and scattered electron and the azimuthal 
angle of the plane define by q and k p . The electron kinematics is contained in 
the kinematical factors 
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The response functions are defined in the standard fashion as 
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In the above expressions, (p( q)) = <f )) denotes the Fourier transform of 
the transition charge density < f\p(r)\i >, while the Fourier transform of the 
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transition three-current density 
'J A Q 
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Current conservation imposes the condition 

q,J^q)=Lup(q)-qJ(q;0) = 0. 
All results presented in this paper are obtained in the Coulomb gauge 

J* 1 = (p, Jx, J y , -P) , 

■ q 

and a relativistic current operator of the form 

J» = Fl (Q 2 )r+^F 2 (Q 2 )a^q v , 



(19) 
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where Fi is the Dirac, F 2 the Pauli form factor and k the anomalous magnetic 
moment. 



3 Relativistic formulation of Glauber theory 



3.1 Dirac-eikonal approximation for potential scattering 



We start our derivations by looking for solutions to the time-independent Dirac 



equation for an ejectile with energy E = ^k^ + M 2 and spin state 
the presence of spherical Lorentz scalar V s {r) and vector potentials V v {r) 
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where we have introduced the notation for the unbound Dirac states. 

ftp ,?Tis 

After some straightforward manipulations, a Schrodinger-like equation for the 
upper component can be obtained 
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where the central and spin orbit potentials V c and V so are defined as 



E V 2 (r) - V 2 (r) 
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Since the lower component w^ + \k p ,m s ) is related to the upper one through 
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the solutions to the Eq. (23) determine the scattering state. In RDWIA ap- 
proaches, a Dirac equation of the type (23) is solved numerically for Dirac opti- 
cal potentials V s (r) and V v (r) derived from global fits to elastic proton scatter- 
ing data [33] . Not only are global parameterizations of Dirac optical potentials 
usually restricted to proton kinetic energies T p < 1 GeV, calculations based on 
exact solutions of the Dirac equation frequently become impractical at higher 
energies. This is particularly the case for approaches relying on partial-wave 
expansions. At higher proton kinetic energies it appears more convenient to 
solve the Dirac equation (23) in the eikonal approximation [20,27,28,34]. In 
intermediate-energy elastic p— 40 Ca scattering (T p ps 500 MeV) the eikonal 
approximation was shown to successfully reproduce the exact Dirac partial- 
wave result. Bianconi and Radici showed that for ejectile momenta exceeding 
1 GeV, the eikonal approximation almost reproduced the 12 C(e, e'p) differen- 
tial cross sections obtained through performing a partial-wave expansion of 
the "exact" scattering wave function [35,36]. 

As shown for example in Refs. [27] and [29], in the eikonal limit the scattering 
wave function takes on the form 
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where the eikonal phase reads (r = (b, z)) 
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In modeling A(e, e'p) processes, the average momentum K occurring in this 
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expression is denned as 



K 



k p + q 
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where q is the three-momentum transfer induced by the virtual photon. The 
scattering wave function from Eq. (26) differs from the plane wave solution in 
two respects. First, the lower component exhibits the dynamical enhancement 
due to the combination of the scalar and vector potentials. Second, the eikonal 
phase e lS ^ accounts for the interactions that the struck nucleon undergoes 
in its way out of the target nucleus. The eikonal approximation is a valid 
one, if the magnitude of the three-momentum transfer | q | is sufficiently 
large in comparison with the projected initial (or, missing) momentum of the 
ejectile (or, q p m — | k p — q |). The eikonal phase of Eq. (27) reflects the 
accumulated effect of all interactions which the ejectile undergoes in its way 
out of the nucleus. All these effects are parametrized in terms of one-body 
optical potentials and the link with the elementary proton-proton and proton- 
neutron scattering processes appears lost. In Glauber theory this link with the 
elementary processes is reestablished. 



3.2 Proton-nucleon scattering 



To start our derivations of a relativistic version of Glauber multiple-scattering, 
we first consider a nucleon-nucleon scattering process and assume that it is 
governed by a local Lorentz and vector potential V s (r) and V v (r). In the eikonal 
approximation, the scattering amplitude corresponding with this process reads 
[27] 
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with a relativistic scattered wave ib? as determined in the eikonal approx- 

T kf,m s / 

imation (Eq. (26)) and the free Dirac solution 
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After some algebraic manipulations one finds [27] 
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where the phase shift function is given by 



X© = IF dz ^ z ) + V °°& z ) P • $ x K)]} . (32) 

A J-oo 

In conventional Glauber theory the phase shift function x(b) is not calculated 
on the basis of knowledge about the radial dependence and magnitude of the 
potentials V c (r) and V so (r), but is directly extracted from proton-proton and 
proton-neutron scattering data. This requires some extra manipulations which 
will be exposed below. 

Using rotational invariance and parity conservation, a pN scattering process 
where at most one polarization of the colliding particles is determined, can be 
analyzed with a scattering amplitude of the form [37] 



V(A) = f; N (A) + f; N (A)a ■ (h x k f ) , (33) 

where f£ N and f£ N are the central and spin-dependent amplitudes, a is the 
spin-operator corresponding with the incident proton, and A = kf — ki is 
the transferred momentum. The small angle elastic scattering of protons with 
k v > 1 GeV is dominated by the central, spin-independent amplitude f£ N 
[37]. Given the diffract ive nature of pN collisions at these energies, the central 
amplitudes are usually parameterized in a functional form of the type 
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The parameters in Eq. (34) can be determined from fitting the results of 
proton-nucleon scattering experiments. A selection of the measured elastic 
(cr el ) and total (a tot ) cross sections for proton-proton and proton-neutron scat- 
tering are shown in Fig. 1. The above form for the scattering amplitude fp N 
corresponds with a differential cross section which reads at forward angles 
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3.3 Relativized Glauber model for A(e, e'p) 



In the relativized Glauber multiple scattering framework, the antisymmetrized 
A-hody wave function in the final state reads 
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Fig. 1. Total and elastic cross sections for proton-proton and proton-neutron scat- 
tering as a function of the proton lab momentum. The data are from Ref. [38]. The 
solid (dashed) curve displays our global fit to the elastic (total) cross section. 
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where ty J A R J^ R is the wave function characterizing the state in which the A — 1 
nucleus is created and A is the antisymmetrization operator. In the above 
expression, the subsequent elastic or "mildly inelastic" collisions of the struck 
proton with the "frozen" spectator nucleons, are implemented through the 
introduction of the following operator 



S (r,f 2 ,r 3 , 
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where r(b,z*j denotes the position of the struck particle and (f 2 ,f 3 , . . . ,r A ) 
those of the (frozen) spectator protons and neutrons in the target. The step 



function 9 (z,- 



guarantees that the struck proton can only interact with 



the spectator protons and neutrons which it finds in its forward propagation 
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path. Further, the profile function for pN scattering is defined as 
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where the last equality can be derived from the expression (34). Two major 
assumptions underly the derivation of the operator of Eq. (37) : the eikonal 
approximation for the continuum wave function of the struck proton and the 
frozen approximation for the positions of the spectators. Indeed, one assumes 
that the ejectile passes through the nucleus in a very short time so that vari- 
ations in the positions of the residual nucleons can be ignored. The operator 
of Eq. (37) represents the accumulated effect of the phase shifts contributed 
by each of the target scatterers as the ejectile progresses through the nucleus. 
The property of so-called phase-shift additivity is a direct consequence of the 
adopted one-dimensional nature of the relative motion, together with the ne- 
glect of three- and more-body forces, recoil effects and longitudinal momentum 
transfer. 

The Dirac-Glauber A(e, e'p) transition amplitude can be written as 
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where for convenience only the spatial coordinates are explicitly written. The 
difficulty of the evaluation of this matrix element stems from the fact that 
the multiple-scattering operator S is a genuine A-body operator. One popular 
approximation in Glauber inspired A(e, e'p) calculations, is expanding the A- 
body operator S (r, r 2 , r 3 , . . . f A ) 
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and truncating it at some order in T. Formally, this expansion bears a strong 
resemblance with the Mayer-Ursell expansion used in modeling the correlation 
effects in the theory of real gases and liquids. In the above expression, the unity 
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operator (first term) reflects the situation whereby the ejectile is not subject 
to scatterings in its way out of the nucleus. The second term, which is linear in 
the profile function, accounts for processes whereby the struck nucleon scatters 
on one single spectator nucleon before turning asymptotically free (single- 
scattering process). Higher-order terms in the expansion refer to processes 
whereby the ejected proton subsequently scatters with two, three, . . ., A-l 
spectator nucleons. Evaluating the different terms in the above expansion 
allows one to estimate the effect of single-, double, triple-, . . . scatterings. 
For heavier nuclei, truncating the expansion at first order in T appears as a 
rather questionable procedure. In Section 4 results obtained with the truncated 
(Eq. (39)) and complete (Eq. (37)) Dirac-Glauber multiple-scattering operator 
will be compared. 

In evaluating the Dirac-Glauber transition amplitude of Eq. (38) we have in- 
troduced a minimal amount of approximations. In line with the assumptions 
underlying the RDWIA approaches we adopt a mean-field approximation for 
the nuclear wave functions. For the sake of brevity of the notations, in the 
forthcoming derivations we consider the case A = 3. A generalization to arbi- 
trary mass number A is rather straightforward. Further, we will assume that 
the nuclear current J M is a one-body operator. As both the initial and final 
wave functions are fully antisymmetrized, one can choose the operator J M 
to act on one particular coordinate and write without any loss of generality 
(f' = (b'(x',y'),z')) 
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with eijk the Levi-Civita symbol, and where we have introduced a frame 
(x', y', z') defined by the following unit vectors 
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The plane (x', z') coincides with what is usually known as the hadron reaction 
plane in A(e, e'p) reactions. 



Adopting a mean-field picture, the initial A-nucleon wave function is of the 
form 
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For spherically symmetric potentials, the solutions Q (r) to a single-particle 
Dirac equation entering this Slater determinant have the form [39] 
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where n denotes the principal, k and m the generalized angular momentum 
quantum numbers. The y± Km are the spin spherical harmonics and determine 
the angular and spin parts of the wave function, 
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The final A-body wave function reads 
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Relative to the target nucleus ground state written in Eq. (42), the wave 
function of Eq. (45) refers to the situation whereby the struck proton resides 
in a state leaving the residual A — 1 nucleus as a hole state in that 

particular single-particle level. 



Assuming that the profile function T does not contain spin-dependent terms, 
one can safely assume that for elastic and mildly inelastic scatterings 
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Inserting this expression in Eq. (40) one obtains 
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This leads to our final result for the Dirac-Glauber A(e, e'p) transition ampli- 
tude 
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where the Dirac-Glauber phase (6, z) is defined in the following fashion 
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In this expression, the product extends over all occupied single-particle states, 
except for the one from which the nucleon is ejected. The RPWIA approxi- 
mation would correspond with putting Q = 1 in the expression for the matrix 
element of Eq. (48). 



The numerical evaluation of the Glauber phase Q (b, z) is rather challenging 
if no additional approximations are introduced. A Monte Carlo integration 
method was suggested in Ref [40]. In our numerical calculations we did not 
introduce any further approximations and found it most appropriate to evalu- 
ate the scattering amplitudes and Glauber phases in the frame defined by the 
unit vectors of Eq. (41). Inserting the expression for the Dirac single-particle 
wave functions <p a of Eq. (43) in the Eq. (49) for the Glauber phase, one gets 

(df' EE dz'b'db'dfa) 
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s$z)= n 



47T/3. 



pN 



db' / cfe'0(,z' - z) 



x 



+ 



— — — ^ Km (i2(6 , z ), a) 



x exp 



A2 



2A 2 



piV J 



2tt 



exp 



A 2 



pN 



(50) 



Standard numerical integration techniques were adopted to evaluate the in- 
tegrals occurring in this equation. It is important to remark that cylindrical 
symmetry about the z' axis makes the above expression to be independent of 
06. As a result the relativistic Glauber depends on only two independent vari- 
ables (b, z). In the above expression (50) each of the frozen spectator nucleons 
is identified by its quantum numbers (n, k, m) and its corresponding Dirac 
wave function nKm (r*, a). In most Glauber-based calculations, an additional 
averaging over the positions of the spectator nucleons is introduced. This pro- 
cedure amounts to replacing in Eq. (50) the characteristic spatial distributions 
of each of the spectator nucleons described by the functions F nK (r) and G m (r) 
by an average density distribution for the target nucleus 



G{b,z) 



^PpN 



b'db'T B (b',z)exp 



A 2 



(b - V ) 

2A 2 



pN 



X 



i-2-ir 

Jo 



V exp 



A 2 



-2sin 



2 <Pb 



pN 



A-l 



(51) 



The function Ts(b' ', z) which was introduced in the above expressions is known 
as the "thickness function" and reads 



1 r+oo 

T B (b',z) = - dz'6(z'-z)p B (r'(b',z')), 
A Jo 



(52) 



where the relativistic radial baryon density p_b(?") is defined in the standard 
fashion 



p B (r) = (n-ro^) = E / (^(r, ^ (^Jr, <?)) 

OLOCC 

(2j + 1 



E 



47rr 2 



\G n Jr) 



\F n Jr) 



(53) 



and the sum over nn extends over all occupied states. 
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Proton momentum [GeV/c] Proton momentum [GeV/c] 



Fig. 2. The Glauber slope parameters (3p p and 0^ n as obtained from Eq. (54) with 
the global fits contained in Fig. 1. The data are from Refs. [41] (proton- neutron) 
and [42] (proton-proton) and are determined from the small-angle t dependence of 
the measured differential cross sections. 

3.4 Glauber parameters 

We refer to A(e, e'p) results obtained with a scattering state of the form of 
Eq. (36) as a relativistic multiple-scattering Glauber approximation (RMSGA) 
calculation. It is worth stressing that in contrast to the RDWIA models, all 
parameters entering the calculation of the scattering states in the RMSGA 
A(e, e'p) model can be directly determined from the elementary proton-proton 
and proton-neutron scattering data. In practice, for a given ejectile's lab mo- 
mentum I k p I the following input is required : the total proton-proton crff* 
and proton-neutron cr*°* cross sections, the corresponding slope parameters 
(Pp P and (3p n ) and the ratios of the real to imaginary part of the scattering 
amplitude (e pp and e pn ). We obtain the numbers a 1 ^, (3p N and € p n through 
interpolation of the data base available from the Particle Data Group [38]. 
The slope parameters (3p p and f^ n may be found by analyzing the shape of 
the differential cross sections assuming that the contribution from the spin- 
dependent terms is negligible. At proton momenta p p < 1 GeV the slope 
parameters found directly from experiment and phase-shift analysis differ sig- 
nificantly due to a large contribution from the spin-dependent scattering am- 
plitude [37]. At higher energies this difference drops quickly indicating that 
spin effects are small in that region. In our calculations, the slope parameters 
are obtained from the ratio of the elastic a pN to the total a P % cross section 
through the following relation 
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10 1 1 10 

Proton momentum [GeV/c] 



Fig. 3. The ratios of the real to imaginary parts of the central amplitude f^ N 
for proton-proton and proton-neutron scattering as a function of the proton lab 
momentum. The curves are our global fits to the data points. The data are from 
the review paper of Ref. [37] . 



In Fig. 2 we compare the slope parameters obtained through this formulae 
with those extracted directly from the t dependence of the differential pN 
cross sections. The curves in Fig. 2 use the above formulae (54) and our global 
fits to cr^, cjp N and € p n shown in Figs. 1 and 3. For proton-proton scattering 
the situation emerges to be very satisfactory. 
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Fig. 4. The radial dependence of the real and imaginary part of the computed 
Glauber phase Q along the direction of the ejected particle (6 = 0°) for proton 
emission from 4 He at various proton kinetic energies T p . Results with the expression 
of Eq. (39) truncated to the first order in T (dashed line) are compared to the full 
result (solid line). 



4 Numerical results for relativistic Glauber phases 



4-1 Single- and multiple- scattering effects 



In many works, in the calculation of the Glauber phases only a limited amount 
of terms in the expansion of Eq. (39) is retained. In some cases, the operator 
S is replaced by the term which is first order in V. This reduces the treat- 
ment of the FSI effects to one-body operators and corresponds with retaining 
the cumulative effect of free passage through the target in addition to sin- 
gle scatterings. We wish to compare Glauber phases obtained with the exact 
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Fig. 5. As in Figure 4 but now for 12 C. 

operator with those that are produced when allowing only single scatterings. 
The Glauber phase, as it was defined in Eq. (49) depends on the two variables 
(b, z) and is independent of <fib- We wish to remind the reader that for conve- 
nience of the numerical integrations the z'-axis is chosen to point along the 
asymptotic direction k p of the ejectile. 

In Figs. 4, 5, 6 results are displayed for the computed real and imaginary part 
of 

Q(r, 9) = Q (b = V r 2 — r 2 cos 2 9 , z = r cos , (55) 

for the target nuclei 4 He, 12 C and 208 Pb and 9 = 0°. Here, 9 denotes the polar 
angle with respect to the axis defined by the asymptotic momentum of the 
ejected particle. We remind that in the absence of final-state interactions the 
real part of the Glauber phase equals one, whereas the imaginary part van- 
ishes identically. As becomes clear from Fig. 4, the approximation of retaining 
only single-rescattering terms, emerges to be a reasonable approximation for 
the 4 He target nucleus. As a consequence, for 4 He the average number of 
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Fig. 6. As in Figure 4 but now for 208 Pb. 

rescatterings can be inferred not to be larger than one. For a given ejectile's 
momentum, the average number of scatterers which it encounters in its way 
out of the nucleus is expected to grow like A 1 ^ 3 . Given that for A = 4 the 
average number of rescatterings is not larger than 1, one can infer that for a 
heavy nucleus like 208 Pb the scattering series of Eq. (39) will receive sizeable 
contributions up to the fourth order in T. This complies with the numbers 
quoted in Table 1 of Ref. [25]. As can be inferred from Figs. 4, 5, 6, single 
rescatterings dominate the real and imaginary part of the Glauber phase at 
the nuclear surface. In the interior of the nucleus, single- and the summed 
higher-order scattering terms come with an opposite sign for all target nuclei 
studied. As a matter of fact, even for a nucleus like 12 C, the truncation to 
single scatterings results in a sizeable overestimation of the FSI effects. The 
real part of the Glauber phase exhibits relatively little T p dependence over 
the range 0.25 - 2.00 GeV covered in the figures. The imaginary part, on the 
other hand, changes sign as one exceeds T v = 0.5 GeV and enters a highly 
inelastic regime. This observed change in the relative sign between the real 
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and imaginary parts of the Glauber phase is governed by the T p dependence 
of the parameters e pp and e pn as is shown in Fig. 3. 

When evaluating the (e, e'p) matrix elements and performing the integrations 
over r and 9, the functions displayed in Figs. 4-6 quantify the effects stemming 
from the FSI along the direction defined by the asymptotic momentum of the 
ejectile. For the Glauber phases, the radial coordinate "r" can be interpreted as 
the distance relative to the center of the target nucleus where the photon hits 
the target nucleus. For a given r, an additional non-trivial integration over the 
polar angles 9 has to be performed. Here, 0° < 9 < 90° (90° < 9 < 180°) refers 
to a situation where the photon hits the nucleon in the forward (backward) 
hemisphere with respect to the direction defined by k p . The dependence of the 
real and imaginary part of the Dirac Glauber phase on the polar angle 9 is 
illustrated in Fig. 7 for emission of 1 GeV protons from 4 He. The 9 = 180° case 
corresponds with a peculiar event whereby the photon couples to the proton 
along the direction defined by —k p . For 9 = 180° and increasing r, the photon 
initially hits the proton at the outskirts of the nucleus and the proton has to 
travel through the whole nucleus before it becomes asymptotically free at the 
opposite side. It speaks for itself that these kinematical situations induce the 
largest FSI effects but cannot be expected to provide large contributions to 
the integrated matrix elements. 



4-2 Relativistic and density effects 

The role of relativity in the description of the FSI can be estimated by ne- 
glecting the small components F nK (r) in the relativistic wave functions for the 
individual scattering nucleons in Eq. (49) and comparing it to the exact result. 
We have performed several of these calculations for a variety of target nuclei. 
One representative result is displayed in Fig. 8. In general, the relativistic 
lower wave function components for the scattering centers (i.e. the nucleons 
residing in the daughter nucleus) are observed to have a minor impact on the 
predictions for both the real and imaginary part of the Glauber phase. In the 
next section, however, it will be shown that inclusion of the lower relativistic 
components is essential for some A(e, e'p) observables. From the results pre- 
sented here, it can be excluded that this could be attributed to a relativistic 
effect in the description of the final-state interactions. 

Besides the neglect of relativistic effects, most Glauber calculations use an 
average nuclear density to describe the spatial distribution for each of the 
frozen nucleons from which the ejectile can scatter. This approximation was 
outlined at the end of Sect. 3.3. One may naively expect that this "averaging" 
at the wave-function level becomes increasingly accurate as the target mass 
number A increases. As a matter of fact, the results displayed in Fig. 9 illus- 
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Fig. 7. The radial dependence of the real and imaginary part of the computed 
Glauber phase Q at various values of 9 and T p =l GeV for proton emission out of 
4 He. 

trate that the quoted approximation is a valid one even for a light nucleus like 
4 He. Inspecting Fig. 9, only in the absorptive part for 4 He and 12 C a minor 
overestimation gets introduced through the averaging procedure. 



5 Numerical results for 4 He(e, e'p) structure functions 



In conventional nuclear-physics models, by which we understand models which 
are based on hadronic degrees-of-freedom, the 4 He nucleus plays a key role. 
Indeed, it is the simplest nuclear system in which all basic characteristics of 
"complex" nuclei exist. Accordingly, it comes at no surprise that the 4 He(e, e'p) 3 H 
reaction plays a pivotal role in investigations into the short-range structure 
of nuclei [43] and medium-modification effects [5]. In Figs. 10 and 11 we dis- 
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Fig. 8. The effect of the lower components in the wave functions for the scattering 
centers on the computed Glauber phase Q at 9 = 0° and T p = 1 GeV for proton 
emission out of 12 C. Results with ignoring the lower components (dashed line) are 
compared to the full result (solid line). 

play some predictions for the separated 4 He(e, e'p) 3 B. response functions in 
kinematics presently accessible at the Thomas Jefferson National Accelera- 
tor Facility As a matter of fact, the kinematics corresponding with Fig. 10 
coincides with a scheduled experiment at this facility [43]. 

Non-relativistic A(e, e'p) calculations typically miss, amongst other things, 
the effect from the coupling between the lower components in the bound and 
scattering states. As a matter of fact, we interpret the contribution from the 
coupling between the lower component in the bound (4> ai ) and scattering 
state (4>k p m s ) to the matrix element of Eq. (48) as a measure for the impact of 
relativistic effects. We denote A(e, e'p) results that are obtained after omitting 
this specific part as "NRMSGA" . We wish to stress that these "NRMSGA" 
calculations still use relativistic kinematics and bound and scattering states 
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Fig. 9. Effect of replacing the wave functions for the individual scattering centers 
by an average density on the computed Glauber phase Q. The results are obtained 
for 9 = 0° and T p = 1 GeV for proton emission out of 4 He, 12 C and 56 Fe. Results 
with the approximated expression from Eq. (50) (dashed line) are compared to the 
exact result from Eq. (49) (solid line). 



obtained by solving a Dirac equation. As can be appreciated from Figs. 10 and 
11, the effect of the coupling between the lower components is marginal for 
the structure functions TZl and TZt, but substantial for the two interference 
functions IZtt and IZtl- These conclusions regarding the role of relativistic 
effects at low missing momenta confirms the major findings of numerous other 
investigations [15,16]. Note that the genuine relativistic effect stemming from 
the coupling between the lower components is very prominent in the (non)- 
filling of the predicted dip in the cross section for p m w 410 MeV. In order 
to illustrate the substantial impact of relativistic effects on some observables, 
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Fig. 10. Separated 4 He(e, e'p) 3 H response functions versus missing momentum in 
constant (q, w)-kinematics under quasi-elastic conditions. The electron kinematics 
was e=4.80 GeV, g=1.50 GeV and cu=0.84 GeV (i.e., Q 2 =1.54 (GeV) 2 ). 

Fig. 12 presents the so-called left-right asymmetry 



corresponding with the kinematics of Figs. 10 and 11. As can be seen the 
genuine relativistic effect due to contribution from the the lower components 
in the bound and scattering state bring about a substantial increase in Alt- 
As a matter of fact, the RMSGA Alt predictions of Fig. 12 are completely in 
line with the optical-potential RDWIA predictions by J. Udias contained in 



The relativistic plane wave impulse approximation (RPWIA) results of Figs. 10 



Alt 



a((f> = 0°) - cr(<f> = 180°) 



vtlKtl 



(56) 




i'l^l + vtR-t + vttTZtt 



Ref. [43]. 
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Fig. 11. Separated 4 He(e,e'p) 3 H structure functions versus missing momentum in 
constant (q, w)-kinematics under quasi-elastic conditions. The electron kinematics 
was e=2.442 GeV, g=1.696 GeV and u=l GeV (i.e., Q 2 =1.88 (GeV) 2 ). 



- 12 are obtained after setting the Glauber phase Q(b,z) equal to one in the 
matrix elements of Eq. (48). This corresponds with a calculation which ignores 
all FSI mechanisms but adopts relativistic kinematics, Dirac bound states, a 
relativistic plane wave for the ejectile and a relativistic current operator. At 
low missing momentum the FSI quenches the cross sections by 20 %, con- 
firming the result of the non-relativistic Glauber calculations of Refs. [20], 
[44] and the calculations of Laget reported in Ref. [45]. Not surprisingly, with 
growing missing momentum the FSI effects become increasingly important. 
For example the dip in the RPWIA response functions at p m ~ 0.5 GeV is 
shifted downwards by about 100 MeV and partially washed out. 
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Fig. 12. The Atl left-right asymmetries as a function of the missing momentum 
corresponding with the kinematics of Figs. 10 and 11. 

6 Summary and conclusions 

A fully unfactorized relativistic formulation of Glauber multiple-scattering 
theory for modeling exclusive A(e, e'p) processes has been presented. For- 
mally, the model bears a strong resemblance with the RDWIA approaches 
which have been developed over the last number of decades. In contrast to 
the RDWIA models, the relativistic Glauber approximation for dealing with 
FSI mechanisms holds stronger links with the elementary proton-proton and 
proton-neutron processes and does not require the (phenomenological) input 
of optical potentials. Our fully unfactorized framework can accommodate all 
relativistic effects which are usually implemented in the RDWIA approaches. 
Like in the RDWIA models, the bound-state wave functions are solutions to 
a Dirac equation with scalar and vector potentials fitted to the ground-state 
nuclear properties, i.e. an approach commonly known as the a — uo model. 

The relative contribution from single- and multiple-scatterings to the inte- 
grated FSI mechanisms, has been investigated. For the target nucleus 4 He, 
it turned out that single-scatterings give already a fair account of the com- 
plete scattering processes, and that double- and triple-scatterings are rare. For 
target nuclei heavier than 4 He, the net effect of double-, triple- and quadruple- 
scatterings is found to affect the real and imaginary parts of the FSI induced 
phase with an opposite sign as compared to the single-scattering term. For 12 C 
and heavier target nuclei, Glauber calculations restricted to single-scatterings 
substantially overestimate the FSI mechanisms. The lower components in the 
relativistic single-particle wave functions, on the other hand, are observed to 
have a negligible impact on the predictions for the distorting and absorptive 
effect of single- and multiple-scattering events which the ejected proton un- 
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dergoes. This observation provides support for the non-relativistic Glauber 
approaches which have been widely adopted in nuclear transparency calcula- 
tions for example. Also the frequently adopted simplification of replacing the 
squared nucleon wave functions for the individual scatterering centers in the 
target nucleus by some averaged nuclear density, was found to lead to results 
which approximate nicely the exact ones. 

The developed RMSGA model has been applied to the exclusive 4 He +e — > 
e'+p+t process. The predicted RMSGA structure functions are found to follow 
similar trends in comparison to predictions from non-relativistic Glauber and 
relativistic optical-potential calculations. Extensions of the presented model 
include the implementation of short-range nucleon-nucleon correlations and 
color transparency effects. Work in this direction is in progress. 
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